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We present a search for the pair production of scalar top quarks, t, using 995 pb _1 of data collected 
in pp collisions with the DO detector at the Fermilab Tevatron Collider at ^/s = 1.96 TeV. Both 
scalar top quarks are assumed to decay into a charm quark and a neutralino (x?), where xl is the 
lightest supersymmetric particle. This leads to a final state with two acoplanar charm jets and 
missing transverse energy. We find the yield of such events to be consistent with the standard 
model expectation, and exclude sets of i and x? masses at the 95% C.L. that substantially extend 
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the domain excluded by previous searches. 
PACS numbers: 14.80.Ly; 12.60.Jv 



Supersymmetry (SUSY) may provide a solution to the 
hierarchy problem if the SUSY particles have masses 
less than 1 TeV, strongly motivating the searches for 
supersymmetric objects at the Fermilab Tevatron Col- 
lider. SUSY predicts the existence of partners with iden- 
tical quantum numbers to all standard model (SM) par- 
ticles except for spin. There exist two spin zero SUSY 
partners of the top quark corresponding to the latter's 
left and right handed states. Several arguments exist 
in favor of a light scalar top quark (t). The t mass 
mj receives negative contributions proportional to the 
top quark Yukawa coupling in the renormalization group 
equations. This makes the t weak eigenstates lighter 
than other squarks Mixing between the left and 

right handed superpartners of the top quark, being pro- 
portional to the top quark mass m ( , leads to a large 
mass splitting between the two physical eigenstates. This 
makes one of the t considerably lighter than the other. 
Additionally, a light t that strongly couples to the Higgs 
boson could also generate a large enough CP violating 
phase to explain the mechanism for electroweak baryoge- 
nesis Q. 

In R-parity conserving models [3j, the lightest super- 
symmetric particle (LSP) is stable, and cosmological con- 
straints indicate that it should be neutral and color- 
less [4J. In the following we assume conservation of R- 
parity and take Xii the lightest of four SUSY particles 
that result from the mixing between the SUSY partners 
of the SM neutral gauge and Higgs bosons, to be the LSP. 

In the search reported in this Letter, we consider the 
range raj < + to^+ and mi < raw + m b + m y° > where 
mfc is the b quark mass, m^o is the \i mass and to~+ is 

the xt mass, with xt being the lighter of two mass eigen- 
states resulting from the mixing of the SUSY partners of 
charged gauge and Higgs bosons. The dominant t decay 
mode in this model is the flavor changing process t — > cxI 
and is assumed to occur with 100% branching fraction. 
The t — > txi decay is kinematically forbidden over the 
t mass range accessible in this search, and the tree level 
four-body decays t — > bff'xi can be neglected @ . 

In pp collisions, t pairs are produced via quark- 
antiquark annihilation and gluon fusion. The t pair 
production cross section (cr—) primarily depends on mj, 
and a weak dependence on other SUSY parameters af- 
fects only the higher-order corrections. At y/s = 1.96 
TeV which is the centre-of-mass energy available at the 
Fermilab Tevatron collider, a— at next-to-leading-order 
(NLO), calculated with PROSPINO 0|, ranges from 15 pb 
to 1 pb for 100 < mj < 160 GeV. These cross sections are 
calculated using CTEQ6.1M parton distribution functions 
(PDFs) 0] and equal renormalization and factorization 



scales /u r f = m^. A theoretical uncertainty of w 20% is 
estimated due to scale and PDF choices. The tt event 
topology consists of two acoplanar charm jets and miss- 
ing transverse energy ($r) from the neutralinos that es- 
cape detection. Searches for t pair production in the jets 
plus missing transverse energy mode have been reported 
by collaborations working at the CERN LEP collider 0] , 
and the CDF [lOU and DO [HE3] collaborations. The 
highest excluded mass to date is < 134 GeV (95% 
C.L.) for m^o = 48 GeV [HI. 

The t search is performed in the data collected with the 
DO detector during Run Ha of the Tevatron and corre- 
sponds to an integrated luminosity of 995 ±61 pb -1 [Hj]. 
A detailed description of the DO detector can be found 
in |l4|. The central tracking system consists of a sili- 
con microstrip tracker and a fiber tracker, both located 
within a 2 T superconducting solenoidal magnet. A 
liquid-argon and uranium calorimeter covers pseudora- 
pidity \r)\ < 4.2, where r\ = — ln[tan(0/2)], and 9 is the 
polar angle with respect to the proton beam direction. 
An outer muon system, covering \rj\ < 2, consists of layers 
of tracking detectors and scintillation counters on both 
sides of 1.8 T iron toroids. 

The data sample collected from April 2003 to Febru- 
ary 2006 with the jets+^r triggers was analyzed for the t 
search. The trigger conditions require the^fr and its sep- 
aration from all jets to be greater than 30 GeV and 25°, 
respectively, where I/It is the transverse energy computed 
only from jets. Jets are reconstructed using an iterative 
midpoint cone algorithm with radius 7?. cono = 0.5 [HI]. 
The data set is reduced to a sample of 1.5 million events 
by requiring at least two jets with px > 15 GeV and 
# T > 40 GeV. 

Signal samples are simulated using pythia 6.323 [Hj 
for mi ranging from 95 GeV to 165 GeV and Xi masses 
from 45 GeV to 90 GeV. The largest expected back- 
grounds for this search are W and Z bosons produced 
in association with jets, denoted as V+jets. The V+iets 
and tt processes are simulated using ALPGEN 2.05 [l7[ 
interfaced with pythia for the generation of initial and 
final state radiation and hadronization. The background 
samples for single top quark and diboson production are 
simulated using COMPHEP [lj| and PYTHIA, respectively. 
The PDF set CTEQ6L1 is used for both signal and back- 
ground samples, and all generated events are subjected 
to full GEANT-based [13] simulation of the detector re- 
sponse. Simulated signal and background events are 
overlaid with recorded unbiased beam crossings to in- 
corporate the effect of multiple interactions that occur in 
a single beam crossing. After reconstruction, simulated 
events are weighted properly to ensure that the instanta- 
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neous luminosity distribution is the same in data and the 
simulated Monte Carlo (MC) samples. A parametriza- 
tion of the trigger efficiency measured from the data is 
applied to simulated MC events in order to fold in trig- 
ger effects. The multijet background, not included in the 
MC samples, is directly estimated from data. 

A large data sample of Z/^*(— * ee) + jets events, cor- 
responding to an integrated luminosity of 1067±65 pb _1 , 
from the same data period as the t search, is used to im- 
prove the prediction of V+jets backgrounds. For this 
study, Z boson candidates are selected using two high 
transverse energy (Et > 15 GeV) clusters that deposit 
more than 90% of their energy in the electromagnetic 
calorimeter, have shower shapes consistent with expecta- 
tions for electrons, are matched with tracks reconstructed 
in the central tracker, and form an invariant mass be- 
tween 65 GeV and 115 GeV. At least two jets with 
Pt > 15 GeV and |f?det| < 2.5 are required, where |?7det| 
is the jet pseudorapidity calculated using the assump- 
tion that the jet originates from the detector center. The 
predicted number of Z/j*(-^ ee)+ > 2 jets events is cal- 
culated using ALPGEN after correcting for differences in 
electron and jet reconstruction efficiencies between data 
and MC and normalizing the MC to the inclusive num- 
ber of Z boson events in data. The ALPGEN prediction 
is corrected in each jet multiplicity bin by a reweight- 
ing function that depends on the transverse momentum 
of the Z boson to obtain better agreement between the 
model and data. The reweighting function is derived by 
fitting the ratio of the transverse momentum distribu- 
tion of Z boson data to that from the MC prediction. 
After reweighting, all other kinematical variables in the 
Z/j*(—> ee)+ > 2 jets sample applicable to the t search 
are well described by MC. 

The multijet background in Z /"/*(— > ee) + jets events 
is estimated from a fit to the dielectron invariant mass 
distribution. The ratio of the number of ee events pro- 
duced by 7* to the number of Z+7* events is determined 
from MC and used to extract the multijet contribution by 
fitting the dielectron invariant mass in data with an ex- 
ponential function for the multijet+7* contribution and 
a Brcit-Wigner convolved with a Gaussian for Z boson 
events. 

For the t search, the predicted SM background from 
V+jets sources is normalized to the number of Z/j*(-^ 
ee) + 2 jets events after subtracting the multijet back- 
ground, ^z{ee)+2' ^ s an example, the normalization 



weight assigned to simulated Z{- 
partons is 



vv) events with n light 



rdata 



T ALP 



W 



Z{vv)+n 
MC 



r Z(ee)+2 u Z{yv)+n e Z{uv)+n 

J srMC 7ALP ~ ' W 

^Z(v9)+n °Z(ee)+21p e Z(ee)+21p 



Here N^P p s +n is the number of simulated Z(—> w)+ 
n light parton jets events; er^l )+n and (?z\L)+2\ P are 
the cross sections predicted by ALPGEN for Z(— > vD)+ n 



and Z/7*(— > ee)+ 2 light parton jets, respectively; and 
£z(vv)+n and ez(ee)+2ip are the corresponding detection 
efficiencies. The factor / = 0.89 ± 0.02 is applied to 
correct for three effects: the absence of 7* contribution 
to Z(vv)-\- jets events, the normalization of MC light jets 
to a data sample that contains all flavors of jets, and the 
difference in the luminosities of the data set used for the 
t search (995 pb _1 ) and the Z/"f*(-^ ee) + 2 jets data set 
(1067 pb- 1 ). 

The normalization weight assigned to simulated W(—> 
£v)+n light partons is 

Ardata _ALP 
, W(lv)+n _ f iV Z(ee)+2 U W(lv)+n e W(tv)+n 



"MC 



where 



= /t7mc Talp 7— a (PT), ( 2 ) 

"W(lv)+n "Z(ee)+21p t Z(ee)+21p 



a{p T ) 



1 d<J w 




r 1 ^ LP i 






' 1 da^V 






<rg LO dp T 


_<?n^ P d P T 



(3) 



is the product of the ratio of the normalized differen- 
tial cross sections for W and Z bosons production at 
NLO 20] and predicted by ALPGEN, respectively. 

The motivation behind using this technique is to lower 
the luminosity times cross section uncertainty (w 6.1%© 
15%) on the predicted number of events towards the 5% 
statistical uncertainty of the Z/j*(-^ ee) + 2 jets normal- 
ization sample. The combined 15% uncertainty on the 
theoretical cross section for various background processes 
is mainly due to the choice of PDF and the renormal- 
ization and factorization scale. The signal and smaller 
backgrounds such as tt, diboson, and single top quark 
production are normalized using the measured absolute 
luminosity. For these processes NLO cross sections were 
computed with MCFM 5.1 



21] 



The search strategy for i involves three steps which 
include the application of the selection criteria on kine- 
matical variables, heavy flavor (HF) tagging and opti- 
mization of the final selection depending on t and Xi 
masses. The data set for_the i search is reduced to a 
sample of 2288 potential tt candidates, by applying the 
15 selection criteria denoted by CI — C15 and summa- 
rized in Table |TJ The main motivation for CI is to reduce 
the multijet background. Requirements C2 to C7 help 
in reducing the VF+jets and multijet backgrounds. The 
condition on the charged particle fraction (CPF) in C8 
requires that at least 85% of the jets' charged particle 
transverse momenta be associated with tracks originat- 
ing from the selected primary vertex in the event. This 
track confirmation requirement removes events with spu- 
rious IpT due to the choice of an incorrect primary vertex. 
C9 — Cll are applied to reject VF+jets background with 
isolated leptons from W boson decay. For an electron to 
be isolated, the energy deposited in the calorimeter in a 
cone of radius 0.4 in rj — cf> around the electron direction 
should not be more than 15% of the energy deposited 
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FIG. 1: Distributions of the asymmetry A = (IpT —^t)/(^t +^t) with the requirement on D = A^ max — A</> m i n inverted (a) 
and applied (c). Distributions of D with the requirement on A inverted (b) and applied (d) for data (points with error bars), 
for SM backgrounds (histogram), and for a signal with m t - = 150 GeV and m^o = 70 GeV (hatched histogram). In all plots 
the signal contribution has been scaled up by five and^?T > 60 GeV is required. The excess in data at A — and D — — 10 
degrees is consistent with the systematic uncertainties on the predicted background. 



TABLE I: Numbers of data events and cumulative signal ef- 
ficiency for mj = 150 and m^o = 70 GeV after each event 
selection. 



Selection 


Events left 


Signal eff. (%) 


Initial selection and trigger 


1.5 x 10 fa 


55.9 


CI: exactly two jets 


464477 


29.5 


C2:$ T > 40 GeV 


440161 


27.5 


C3: A^QetjJeta) < 165° 


278505 


26.5 


C4: jet-1 p T > 40 GeV 


216382 


24.7 


C5: jet-1 T7 do t < 1.5 


113591 


24.6 


C6: jet-2 p T > 20 GeV 


80987 


22.0 


C7: jet-2 |r, dct | < 1.5 


62910 


20.1 


C8: jet-1 jet-2 CPF > 0.85 


49140 


19.8 


C9: isolated track veto 


23832 


13.4 


CIO: isolated electron veto 


23194 


13.3 


Cll: isolated muon veto 


23081 


13.3 


C12: A</> max - A<?Wn < 120° 


9753 


12.6 


C13: A > — 0.05 


3733 


12.0 


C14: A<Xjet, $t) > 50° 


3375 


11.6 


C15: # T > 60 GeV 


2288 


10.0 



in the electromagnetic layers inside a cone of radius 0.2. 
A muon is declared isolated if the sum of the energies 
of all tracks other than the muon in a cone of radius 
0.5 is less than 2.5 GeV, and the calorimeter energy de- 
posited in a hollow cone with inner and outer radii 0.1 
and 0.4 around the muon direction is less than 5 GeV. A 
track with px > 5 GeV is considered isolated if no other 



track with px > 1.5 GeV is found in a hollow cone of 
inner and outer radii 0.1 and 0.4 around the track con- 
sidered. This condition also helps suppress backgrounds 
with t leptons where the r decays hadronically. Remain- 
ing instrumental background is removed using a quan- 
tity defined by the angular separation between all jets 
and the Ifix of the event, D = A0 max — A0 m j n , where 
A<?i max (A(/> m i n ) is the largest (smallest) azimuthal sepa- 
ration between a jet and-^T; and an asymmetry variable 
defined as A = (ffir -$t) I ($t +$t)- The requirements 
applied on these variables are given by C12 and C13. 
Figure Q] shows that both of these variables are very effec- 
tive in eliminating multijet background which dominates 
in data for large D and negative A. 

The 2288 events selected in data can be compared to 
the 2199 ± 18^321 events predicted from the simulation 
normalized to Z/j*(—> ee) + 2 jets events or 2292±19l532 
events predicted using absolute luminosity normaliza- 
tion, with the first quoted uncertainty due to finite MC 
statistics and the second due to systematic effects de- 
scribed in more detail below. The small remaining mul- 
tijet background in the t search analysis is estimated af- 
ter applying all analysis conditions shown in Table U ex- 
cept that onlpT- The background subtracted $r distri- 
bution is fitted in the control region (40 < $t < 60 GeV) 
with exponential and power law functions, and the esti- 
mated contribution is extrapolated into the signal region 
(ffix > 60 GeV) . The average of the two results is taken 
as the multijet background estimate, while the difference 
between the two fit results is taken as the systematic 
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FIG. 2: Distributions of S = A^ max + A0 m j n for data (points 
with error bars), SM background (histogram), and a signal 
with rt%i — 150 GeV and m^o = 70 GeV (hatched histogram) 
after requiring HF tagging but before optimization. 



uncertainty. This amounts to 14.4 ± 10.7 (stat) ± 5.1 
(sys) events contributed by multijet background before 
HF tagging and optimization of selection cuts. 

After selecting candidate events on the basis of topol- 
ogy, HF tagging is used to identify charm jets in the final 
state. A neural network (NN) tagging tool 22j that com- 
bines information from three different DO HF taggers to 
maximize the b quark tagging efficiency (sa 73%) is used 
for this purpose. 

The first tagger converts information from the impact 
parameter of the tracks identified in a jet into a proba- 
bility that all tracks originate from the primary vertex, 
where the impact parameter is the distance of closest 
approach to the interaction point in a plane perpendic- 
ular to the beam axis. The second tagger identifies the 
presence of vertices that are significantly displaced from 
the primary vertex and associated with a jet. The third 
tagger makes use of the number of tracks with large im- 
pact parameter significance, where the significance is de- 
fined as the ratio of the impact parameter to its uncer- 
tainty. The result of the combination is a NN output. 
A requirement on the NN output is made that preserves 
high efficiency for detection of charm jets (« 30%) with a 
ss 6% probability for a light parton jet to be mistakenly 
tagged. The efficiency for c jet tagging is obtained by 
scaling the b jet tagging efficiency measured in the data 
by the c-tagging-to-6-tagging efficiency ratio computed in 
the MC. 

At the final stage of the analysis, additional selection 
criteria on three kinematic variables; $t, S = A</> max + 
A(/> m ; n , and Ht, with Ht defined as the scalar sum of the 
Pt of all jets, are optimized by maximizing the expected 
lower limit on the neutralino mass for a given m t ~. The 
variable S after requiring at least one jet in the event to 
be HF tagged is shown in Fig. [2] 

Minimum values of Ht are varied from 60 GeV to 
140 GeV in steps of 20 GeV, while those for^r are var- 



TABLE II: Optimized values of selections, numbers of ob- 
served data and predicted background events. A requirement 
oilpT > 70 GeV was chosen in all cases. The values of m,- 



and Ht are 


in GeV while those for S are in 


degrees. 


'"? 


H T 


S 


Observed 


Predicted 


95 - 130 
135 - 145 
150 - 160 


> 100 

> 140 

> 140 


< 260 

< 300 

< 320 


83 
57 
66 


85.3±1.8±^ 
59.0 ± 1.6±|;| 
66.6±l.lt? 6 



TABLE III: For three t and Xi ma ss combinations: signal 
efficiencies and the numbers of signal events expected. The 
first errors are statistical and second systematic. The nominal 
(NLO) signal cross section and upper limit at the 95% C.L. 
are also shown. 





Efficiency 


Expected Signal 


& n om 


0"95 


GeV 


(%) 


Events 


P b 


pb 


(130,55) 


1.5 


51.9 ± 2.7t 7 7 i 


3.44 


2.41 


(140,80) 


0.9 


19.6 ± o.stli 


2.24 


2.87 


(150, 70) 


2.1 


30.8 ± 1.2+" 


1.49 


1.42 



ied from 60 GeV to 100 GeV in steps of 10 GeV. Events 
having the values of these quantities above the minima 
are kept. Maximum values of S are tested between 240° 
and 320° in steps of 20°, and events having S below the 
minimum are retained. For each set of requirements, the 
expected value of the signal confidence level (CL S ) Q 
under the hypothesis that only background is present is 
evaluated using all t and Xi mass combinations, taking 
into account systematic uncertainties. The set of criteria 
that return (CL S ) = 0.05 for the highest neutralino mass 
corresponding to a given are chosen to be the optimal 
ones. 

The optimized values of the selections for different rrif 
are given in Table [II] along with the number of events ob- 
served in data and expected SM background. In all cases 
a requirement of|?T > 70 GeV is imposed. No contam- 
ination remains from multijet background at this point 
in the analysis; it is therefore neglected while setting the 
limit. Efficiencies for three signal mass points along with 
the expected numbers of events are shown in Table IIIII 
The distribution of Ht after optimization but with the 
constraint on Ht removed is shown in Fig. (3) The final 
distribution of$r is shown in Fig. 2] The detailed SM 
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FIG. 3: Distributions of Ht after applying optimized require- 
ments on IpT and S for data (points with error bars), SM 
background (histogram), and a signal with = 150 GeV 
and m-o = 70 GeV (hatched histogram). 
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FIG. 4: Final distributions oilpT for data (points with error 
bars), SM background (histogram), and a signal with mj = 
150 GeV and m^o = 70 GeV (hatched histogram). 



background composition is given in Tabic HVl 

Systematic uncertainties are evaluated for each t and 
X\ mass combination for the optimized set of require- 
ments. Sources of systematic uncertainty include jet en- 
ergy scale, jet energy resolution, jet identification and re- 
construction, the jet multiplicity requirement, trigger ef- 
ficiency, data to MC scale factors, normalization of back- 
ground, HF tagging, luminosity determination, choice of 
PDF, and W boson pt reweighting. The effect of the 
jet multiplicity requirement on the background is stud- 
ied using Z/^*(— ► ee) + jets events. The spectrum of 
transverse momentum of the third jet in data events with 
three or more jets is observed to be very well described by 
the simulation generated with ALPGEN. The » 1% sta- 
tistical uncertainty of the lowest px bin, where the bulk 
of the events are, is taken as a systematic uncertainty in- 
troduced by the jet multiplicity requirement. To study 
the effect of the same requirement on the t signal, where 



TABLE IV: Numbers of predicted background events from 
different SM sources for a selection optimized for m t - > 
150 GeV. The uncertainties are due to the limited MC statis- 
tics. 



SM process 


Number of events 


W(-» Iv) + jets 


20.0 ± 0.7 


Z{— * vv) + jets 


15.8 ± 0.5 


tv) + HF (66, cc) 


12.6 ± 0.5 


Z(-» vv) + HF (66, cc) 


11.6 ± 0.4 


ti and single top 


3.7 ± 0.1 


WW, wz, zz 


2.7 ± 0.1 


Z(-> U) (e,/j,,T) + jets 


0.1 ± 0.01 


Z(-> 11) \e,n,r) + HF (66, cc) 


0.1 ± 0.01 


Total 


66.6 ± 1.1 



a third jet enters an event primarily through initial or 
final state radiation, the pt spectrum of the leading jet 
in simulated Z/j*(— > ee) events generated with PYTHIA 
is examined. Comparison between data and simulation 
shows a slight excess in data in the low px bin; this dis- 
crepancy is used to estimate a systematic uncertainty of 
±1.5% on the signal acceptance attributable to the jet 
multiplicity requirement. The uncertainty on the signal 
acceptance and background estimation due to the PDF 
choice was determined using the CTEQ6.1M PDF set. 

The combined 10% uncertainty on the background nor- 
malization includes: 5% uncertainty from Z/j*(— ► ee) + 
jets statistics assigned to all V+jets samples; 50% uncer- 
tainty on the NLO cross section assigned to the V + HF 
background; 6.1% luminosity uncertainty assigned to ti, 
diboson, and single top quark background; and 8%, 6% 
and 15% uncertainties on NLO cross sections for ti, dibo- 
son, and single top quark production, respectively. The 
uncertainty on the background estimation due to the W 
boson pt reweighting is estimated using two different 
methods to estimate the VF+jets background. In the 
first method, the W+jets background is estimated using 
the expression given in Eq. O In the second method, 
the same reweighting function as applied to the Z boson 
was used to reweight the W boson px which is equivalent 
to setting a{px) = 1 in Eq. [2] Detailed estimates of all 
systematic uncertainties are given in Table IVl 

Using the assumption that t decays into a charm quark 
and a neutralino with 100% branching fraction and the 
nominal t pair production cross section, the largest 
excluded by this analysis is 155 GeV, for a neutralino 
mass of 70 GeV at the 95% C.L. With the theoretical 
uncertainty on the t pair production cross section taken 
into account, the largest limit on is 150 GeV, for 
m^o = 65 GeV. These results are shown in Fig. [5] 

In summary, DO has searched for scalar top quarks 
in jets plus missing transverse energy final states using 
1 fb _1 of data. No evidence for t production has been 
found. This analysis substantially extends the excluded 
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TABLE V: Breakdown of systematic uncertainties on the SM 
background and for a signal point with m t - = 150 GeV and 
m.o = 70 GeV. 



Source 


SM 

background 


Signal 


Jet energy 


+ 1.7% 


±2% 


scale 


-2.5% 


-4% 


Jet resolution 


±1% 


±1% 


Jet reconstruction 


±0.8% 


±0.1% 


and identification 






r T"Vl O'crpv 

±11 6&ei 


±6% 


±6% 


Scale factor 


±5% 


±5% 


Normalization 


±10% 




Luminosity 




±6.1% 


HF tagging 


±4.1% 


±3.5% 


PDF choice 


±4% 


±8.7% 






-5.5% 


Two jet cut 


±0.9% 


±1.5% 


W boson pt 


±3% 




reweighting 
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FIG. 5: Region in the t- Xi mass plane excluded at the 95% 
C.L. by the present search. The observed (expected) exclusion 
contour is shown as the green solid (dashed) line. The yellow 
band represents the theoretical uncertainty on the scalar top 
quark pair production cross section due to PDF and renormal- 
ization and factorization scale choice. Results from previous 
searches M, LLCj, Q2l are also shown. 



region of the t — Xi mass plane over the searches carried 
out previously. 
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